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Abstract 

Recently Richard Stanley initiated a study of the distribution of the length 
asn{w) of the longest alternating subsequence in a random permutation w from 
the symmetric group 5„. Among other things he found an explicit formula for 
the generating function (on n and k) for Pr (asn{w) < k) and conjectured that the 
distribution, suitably centered and normalized, tended to a Gaussian with variance 
8/45. In this note we present a proof of the conjecture based on the generating 
function. 

If w = {wi ■ ■ ■ Wn) is a permutation in the symmetric group iS„ then an increasing 
subsequence of length A; is a subsequence Wi-^ ■ ■ ■ Wi^, satisfying 

Wi^ < Wi^ - ■ ■ < Wi^. 

The random variable is„(w), the length of the longest increasing sequence in a random 
permuation w from iS„, has been much studied. Its mean was first determined asymp- 
totically by Logan-Shepp |2] and Vershik-Kerov j3], proving a conjecture of Ulam, and 
the limiting distribution was determined in the celebrated work of Beik-Deift- Johansson 

Recently Stanley jH] initiated a study of the distribution of the length of the longest 
alternating subsequences of w, a subsequence such that 

Wi^ > < Wi,^ > ■ ■ ■ Wi^ . 

If asn{w) denotes the length of the longest alternating sequence in a random permutation 
w from Sn denote by Pn{k) the distribution function for as„(w), 

Pn{k) = : SiSn{w) < k}. 
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Stanley found an explicit formula for the generating function 

B{x,y) = J2 Vn{k)y^x''. 

k,n>0 

It is given by 

, l+p + 2?/e^" + (l-p)e2''- 
B{x,y) - 



1 + p - ?/2 + (1 _ p _ y2) g2px ' 

where p = \/l — y'^. Using certain consequences of this formula he made the conjecture 
that the limit 

Kit) = lim pJ2n/3 + y/nt). 

n~*oo 

exists. We shall show here that this is so, and that K{t) is a Gaussian with variance 
8/45: 



K{t) = ^ / e-' ds. 

a/tT J-oo 

Two other proofs of the conjecture, very different from this one, are indicated in 
We use the integral representation 

Pn{k) = ttAt? / / B{x,y)x-''-^y-''-Uxdy, 



{2TTi) 

where the contours of integration are little curves around zero. We begin by making the 
substitution 

y = sech u, 

where u runs along a vertical line from a + ire to a — in, with a a large positive real 
number. Then y runs counterclockwise around a little curve surrounding zero. 

With this substitution 

p = Vl — sech = tanh u 

(when y is small p is close to 1, which is why p is given by this square root rather than 
its negative) and we find that 

r r fi^^ 2, 6^ 6^ tanh u _j_ g2 x tanh u 

^-^^^ = J2^-I I e^--e^-tanh. « dx du, 

where now u runs from a — ni to a + ni. 

We perform the x integration first. The poles of the integrand are at x = and 

X = {u + jni) cothu, — oo < j < oo. 

The residue of the quotient in the integrand at the jth point is computed to be zero 
when j is odd (so there is actually no pole there) and — 2cothM when j is even. We 
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evaluate the x-integral by integrating around a sequence of expanding contours passing 
half-way between the poles^ and find that 



If °° 
p^(k) = — / cosh'^ti tanh"ti V {u + 2j7riy~^du. (1) 

j=-oo 

(Since j was even we replaced it by 2j.) 

Because of the periodicity of the integrand, we may take any real a and as contour 
of integration any curve from a — ni to a + ni as long as it stays in the region |Imti| < tt 
and to the right of the poles at zero and ±i7r/2. 

Here is how we will proceed. Because u lies in the region |ImM| < tt it seems likely 
that the term of the sum corresponding to j = will dominate the rest. Suppose this is 
so, and that k will be of the order n, say equal to an + o{n) for some positive constant 
a. Then the main part of the resulting integrand will be 

{coshur (^)^ 

which wc write as e"'^'^"\ The method of steepest descent, or saddle point method, tells 
us that we should try to take as our contour of integration one on which Im a{u) is 
constant and on which Rea{u) achieves its maximum at a point Uq. If Uq is not an end- 
point of the curve then a'{uo) — 0. This uq is a saddle point for the function \e"'^'^^\ and 
at all points of the curve its direction away from the saddle point is that of most rapid 
decrease of the function — whence the names for the method. The main contribution to 
the integral will come from the immediate neighborhood of the single saddle point, and 
so one only has to use an expansion of the integrand near this point to determine the 
asymptotics. In practice one first locates a saddle point, follows a steepest descent curve 
until it ends (at a zero of the integrand, at another saddle point, or the boundary of the 
region which could include oo), and then sees if the original contour may be deformed 
to it. 

In our case the saddle points are the zeros of 

f sinhii 1 1 

coshii sinhii coshii u 

We shall show that there are three saddle points in the region \lmu\ < tt, and that they 
coincide when a — 2/3. 

^The quotient in the integrand equals 1 plus a constant times (1 + e''"^ *'*"'^")/ sinli(M — x tanhw). 
There is an expanding sequence of contours on which this is uniformly bounded. A linear change of 
variable shows that it is enough to show this for the function (1 + e^)/ sinh2;. One can take the jth 
contour to be the square with vertices ±{j + 1/2) tt ± + 1/2) ni. 
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Let us consider Aargcr'(u) over the boundary of the region, with httle semi- circular 
indentations of diameter 6 below the pole at irr and above the pole at — ztt. Using the 
fact that cr'('u) is odd we need consider only the upper part. We also use the fact that 
when Im u = it the first two summands in the expression for (j'{u) are real while —1/u 
has positive imaginary part. As u goes from +00 + in to — 00 + ivr along the upper 
boundary with indentation cr'{u) goes from a to 6^^ + 0(1) staying in the upper half- 
plane, then counterclockwise around a large semi-circle (more precisely, within 0(1) of 
one) to —6~^ + 0(1), then to —a while again staying in the upper half-plane. Clearly, 
then, A argcr'(M) over this upper boundary is vr, and so A a.Tga'{u) over the full boundary 
is 2n. This means that in our region the number of zeros of ^'(u), counting multiplicity, 
is one more than the number of poles. There are two poles, at ±m/2, so there are three 
zeros. 

We compute that in the neighborhood of m = 

aiu) = a - ^u' + 0(««)) - + ^u' + Oiu% 

When a = 2/3 this has a zero of order four at m = 0, so cr'{u) has a zero of order three 
and therefore this is the only saddle point. So we do take a = 2/3, and near u = 

The next step is to determine the steepest descent curve, or at least to describe it 
in sufficient detail. It is a curve (or curves) emanating from the saddle point on which 
lma{u) is constant (in this case zero) and on which lie a{u) decreases as we move away 
from the saddle point. We see that there are four such curves, emanating from in the 
directions arg u = ±7r/4 in the right half-plane and the directions arg u = ±37r/4 in the 
left half-plane. By the symmetries of cr{u) (namely that it is even and takes conjugate 
values at conjugate values of u) these curves are reflections of each other in the real and 
imaginary axes. On all of them cr{u) is real and decreases as we move away from the 
saddle point. None of the four curves, when extended, can intersect in the interior of 
the region in question because if they did it would have to be at a saddle point other 
than u = 0, and there are none. It follows from this, and the symmetry of the four 
curves, that the ones emanating from zero in the directions arg u = ±7r/4 stay in the 
right half-plane and meet the upper boundary of our region at points a ± ivr with a > 0. 
The two together, from a — iirtoOtoa + iTT, constitute the curve we use as our contour. 
We denote it by F. Its important properties are its location (to the right of the poles at 
±i7r/2 and therefore giving the same integral as the original contour) and the fact that 
on it (t{u) is real and decreases from zero as we move away from the saddle point. 

Let us now take k = 2n/3 + \pnt in ((T)). The integrand has a pole at m = 0, coming 
from the term of the sum with j = 0. Since was to the left of the original contour of 
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integration, we replace the term 1/u by — e), then deform the contour to F, and 
then take the hmit as e 0. This is indicated by replacing the factor 1/u by — 0). 
(The integral arising from this term can also be expressed in terms of a principal value 
integral.) Thus (0) becomes in present notation 

du 



1 c (ill 

pj2n/3 + v^t) = — / e"'^^") (coshM)v^* 

m Jr u — 

(coshu)v^* J2 



du 



— / e 

TTt Jr 



^^\u + 2j'Ki) u + 2jTTi 



We shall use the notations c and C to denote constants, small and large respectively, 
which will vary with each use. 

We first show that the integral involving the sum is exponentially small. The curve 
r lies in some boundeed subregion of |ImM| < vr. If \j\ is sufficiently large then we shall 
have 



u 



u + 2j7ii 



< 



1 



uniformly on F, and the sum over these j will be 0(2 "). Each of the finitely many 
other summands has absolute value at most 



u 



u ± 27ii 

the sign being that of j, because u ± 27ri and u + 2j7rz have the same real part but the 
absolute value of the latter is larger than that of the former when |j| > 1. The other 
factor in the integrand, {u + 2j7ii)~^, is bounded. It follows that the absolute value of 
the integral is at most a constant times 



^ e''"(") |coshM|v^* 


u 


n 

+ 


u 


n N 

+ 2-"^] \du\ 


U + 27Ti 


u — 27Ti 







In a sufficiently small neighborhood of zero in T the expression is parentheses will be 
0(e^^"). The cosh factor is 0(6*^^) and the first factor is at most 1. Thus the product 
is 0{e~^"'). Outside this neighborhood e"'^*^"-' = 0(e~'^"') because a{u) decreases from 
zero away from the saddle point. The expression in parentheses is 0(1) and the cosh 
factor is as before. Hence the integrand is 0(e~'^") on F and therefore so is the integral. 

We have shown that the integral involving the sum is exponentially small. Let us 
now consider the other integral. 

Take any 6 > and denote by Fi the part of F on which \u\ < and by F2 
the remainder of F. Then e"°"*^"^ = 0{e~'^^ ) on F2, for this is the estimate when 
\u\ = and e"'^^'"^ is less than this on F2. (Notice that arg-u^ is close to ivr when 
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|m| = n Including the other factors in the integrand, we see that on r2 the integrand 
is 0(e~'^"^ ^^'^'-''^^^^ n^). If we choose 5 < 1/8, and we do, then this bound is just 
Q^Q-cn y Hence the integral over r2 is exponentially small. 

Finally we come to the part that gives the asymptotics, the integral over Fi. We use 
the expansion 

e""(") (coshM)v^* = exp [n {u^Ab + 0{u^)) + {u^ /2 + 0{u^))] 

valid near m = 0. We want to replace Fi by F which consists of the portion of the rays 
argu = ±7r/4 from zero to the points where \u\ = n^^ . The difference between the two 
integrals is the integral over the line segment joining the ends of Fi and F, on which the 
integrand is exponentially small. So, with exponentially small error, we may make this 
replacement. The contribution of the integral over F may be written 

— / exp in iu^lAh + Oiu'')) + + Oiu^))]^^. 

Tii Jr ^ J M — 

The change of variable u = rT^l^v^l'^ replaces this by 

-L / f/0/._ ^/ i/oQ\\ ./ l/2 2v 



— r exp f(t;V45 + 0(72-1/2 ^=^))^w/2 + 0(n-i/2 



v-0 



' exp{{vyA5 + 0{n-^^v^))+t{v/2 + 0{n-^^v))}- 



2711 J-in^/2-2S L J f — 

The exponential factor is uniformly bounded by an integrable function e-^l^'l^+'^l"! 
so^ we can take the limit under the integral and find that the above has limit 

1 r'"^ ^vVi5+tv/2 dv 



2m J-ioo V — 

This is equal to 

—= / e " ds. 

J-oo 
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^This is not quite enough. An integral fn{v) dv/{v — 0) is equal to 

TTi /„(0) + PV / fn{v) dv/v = ni /„(0) + / (/„(«) - /„(-«)) dv/v. 



To show convergence to the limit f{v) dv/{v — 0) it is enough to have pointwise convergene of /„ 
to /, dominated convergence outside a neighborhood of and, for example, uniform boundedness of 
the derivatives in this neighborhood. This is seen to hold in our case. 
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